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Abstract. We calculate the exact Casimir interaction energy between two perfectly 
conducting, very long, eccentric cylindrical shells using a mode summation technique. 
Several limiting cases of the exact formula for the Casimir energy corresponding to this 
configuration are studied both analytically and numerically. These include concentric 
cylinders, cylinder-plane, and eccentric cylinders, for small and large separations 
between the surfaces. For small separations we recover the proximity approximation, 
while for large separations we find a weak logarithmic decay of the Casimir interaction 
energy, typical of cylindrical geometries. 
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1. Introduction 

Almost 60 years ago pQ, Casimir discovered an interesting macroscopic consequence 
of the zero point fluctuations of the electromagnetic field: an attractive force between 
uncharged parallel conducting plates. Since then, the dependence of the Casimir force 
with the geometry of the conducting surfaces has been the subject of several works 
[2|. For many years, the only practical way to compute the Casimir energy for non 
planar configurations was the so called proximity force approximation (PFA) [Hj. This 
approximation is valid for surfaces whose separation is much smaller than typical local 
curvatures. Due to the high precision experiments performed since 1997 [3], there 
has been a renewed interest in the geometry dependence of the Casimir force, and 
in particular in the calculations of the corrections to the proximity approximation. 

In the last years there have been a number of attempts to compute the 
Casimir forces beyond the PFA, using for example semiclassical j3J E] and optical [7] 
approximations, and numerical path- integral methods [Hj. Large deviations from PFA 
for corrugated plates have been reported j^j, and in recent months the Casimir energy 
has been computed exactly for several configurations of experimental interest, as the 
case of a sphere in front of a plane, and a cylinder in front of a plane ^HJEIEIEI!- As 
first suggested in [Tl], the latter configuration is intermediate between the sphere-plane 
and the plane-plane geometries, and may shed light on the longstanding controversy 
about thermal corrections to the Casimir force. There is an ongoing experiment to 
measure precisely the Casimir force for this geometry [TH] . 

The configuration of two eccentric cylinders is of experimental relevance too 
|141 ITol IT7] . Although parallelism is as difficult as for the plane-plane configuration, 
the fact that the concentric configuration is an unstable equilibrium position opens 
the possibility of measuring the derivative of the force using null experiments. Up to 
now, the Casimir interaction energy between two cylindrical shells has been computed 
semiclassically and exactly in the concentric case jHl EE], and using the proximity 
approximation in the eccentric situation ^3J QS]. In principle, one could consider 
experimental configurations in which a very thin metallic wire is placed inside a larger 
hollow cylinder. In this case, a more accurate determination of the Casimir force is 
needed. The aim of this paper is to describe in detail the derivation of the exact Casimir 
interaction energy for eccentric cylinders, initially reported by us in J7| , and to compute 
analytically different limiting cases of relevance for Casimir force measurements in this 
configuration. To this end we will use the mode summation technique combined with 
the argument theorem in order to write the Casimir energy as a contour integral in the 
complex plane [TH] . 

The paper is organized as follows. In Section 2 we derive an expression for the 
Casimir interaction energy for any configuration invariant under translations in one of 
the spatial dimensions. When properly subtracted, this expression reduces to an integral 
over the imaginary axis, and is similar to expressions for the Casimir energy derived 
using path integrals or scattering methods. In Section 3 we derive the exact formula for 
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the interaction energy between eccentric cylinders and we analyze some particular cases 
of the exact formula. We first show the known results for the concentric case obtained 
from the exact formulation, and that it is possible to derive the interaction energy for 
the cylinder-plane configuration in the appropriate limit. In Section 4 we consider the 
exact formula in the limit of quasi concentric cylinders of arbitrary radii. We discuss two 
opposite limits of this exact formula: large and small separations between the metallic 
surfaces. In the first limit, we find that the Casimir energy between a thin wire contained 
in a hollow cylinder has a weak logarithmic decay as the ratio between the outer and 
inner radii becomes very large. In the second limit, we recover previous results obtained 
using PFA for quasi concentric cylinders. Finally, Section 5 contains the conclusions of 
our work. 



2. Casimir energy as a contour integral 

The Casimir energy for a system of conducting shells can be written as 

E c=7jY,( w p-™p)> (!) 
z p 

where w p are the eigenfrequencies of the electromagnetic field satisfying perfect 
conductor boundary conditions on the surfaces of the conductors, and w p are those 
corresponding to the reference vacuum (conductors at infinite separation). Throughout 
this paper we use units h = c = 1. The subindex p denotes the set of quantum 
numbers associated to each eigenfrequency. Introducing a cutoff for high frequency 
modes E c (o~) = \ Yjp{ e ~ aWpw p ~ e~ aWp w p ), the Casimir energy E c is the limit of E c (o~) 
as a — > 0. For simplicity we choose an exponential cutoff, although the explicit form is 
not relevant. 

Let us consider a general geometry with translational invariance along the z— axis 
(as for example very long and parallel waveguides of arbitrary sections). The transverse 
electric (TE) and transverse magnetic (TM) modes can be described in terms of two 
scalar fields with adequate boundary conditions. In cylindrical coordinates, the modes 
of each scalar field will be of the form h nj k z (t, r, 9, z) = e^~ lWn ^ t+lkzZ ^ R n (r, 9), where the 
eigenfrequencies are w U)kz = \Jk 2 z + A^, and \ n are the eigenvalues of the two dimensional 
Laplacian 

(£ + ;| 4£ + *)*■<••■ s > = a < 2 > 

The set of quantum numbers p is given by (n, k z ). For very long cylinders of length L 
we can replace the sum over k z by an integral. The result is 

E ^ a) = i£LIf£ ( V / ^ + ^ e " CTv/ ^ f " V^T^e-v^) .( 3 ) 



From the argument theorem it follows that 

— [ d\ \ e-° x 4r 
2m Jc d\ 



h L dX x e ~ aX ^ ln /( A ) = E A < e ~ aXi > ( 4 ) 
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Figure 1. Geometrical configuration studied in this paper. Two perfectly conducting 
eccentric cylinders of radii a < b, length L, and eccentricity e interact via the Casimir 
force. The figure on the right shows the polar coordinates (r, 9) and (p, ip) of any 
point P between the eccentric cylinders used for the determination of the classical 
eigenvalues for this configuration. 

where /(A) is an analytic function in the complex A plane within the closed contour C, 
with simple zeros at A 1; A 2 , • • • within C. We use this result to replace the sum over n 
in Eq.fjHJ) by a contour integral 




where F is a function that vanishes at X n for all n (and F^ vanishes at A n ) . 

In the rest of this section we will consider the particular configuration of two 
eccentric cylinders with circular sections of radii a and b, respectively. We will denote 
the eccentricity of the configuration by e (see Fig. 1). The geometrical dimensionless 
parameters a = b/a and 5 = e/a fully characterize the eccentric cylinder configuration. 
It is worth emphasizing that the results of this section can be trivially extended to more 
general configurations, as long as they are translationally invariant along one spatial 
dimension. It is convenient to compute the difference between the energy of the system 
of two eccentric cylinders and the energy of two isolated cylindrical shells of radii a and 
b, 

E 12 (cr) = E c (a) - E 1 (a, a) - E^cr, b) , (6) 

where 




Here -F lcy i(a) is a function that vanishes at the eigenfrequencies of an isolated cylindrical 
shell of radius a. Therefore 




where 

F Fi^(oa) 2 
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Figure 2. Contour for the integration in the complex plane. 



To proceed we must choose a contour for the integration in the complex plane. 
In order to compute E c (a), E 1 (a, a), and E 1 (a,b) separately, an adequate contour is a 
circular segment Cr and two straight line segments forming an angle and it — with 
respect to the imaginary axis (see Fig. 2). The nonzero angle is needed to show that 
the contribution of Cr vanishes in the limit r —>■ oo when a > 0. For the rest of the 
contour, the divergences in E c (a) are cancelled out by those of Ei(a, a) and Ei(a, b), as 
in the case of concentric cylinders [BJ. Therefore, in order to compute Ei 2 (a) we can set 
= and o" = 0, and the contour integral reduces to an integral on the imaginary axis. 



As we will see, M(iy) is a real function - hence, the integral over y in Eq. (fTU|) is 
restricted to y > k z . After some straightforward steps one can re- write this equation as 



As we have already mentioned, a similar expression can be derived for conductors of 
arbitrary shape, as long as there is translational invariance along the z-axis. It is worth 
noting that the structure of this expression is similar to the ones derived recently for the 
cylinder-plane geometry using path integrals fTJ ^] , and for the sphere-plane geometry 
using the Krein formula 10 . 

3. The exact formula 

In this section we derive the exact formula for the Casimir energy between eccentric 
cylinders. We proceed in two steps: we first find the function F with zeroes at the 
eigenfrequencies for the geometric configuration. Then we obtain an explicit expression 
for the function M, which involves a definition of the Casimir energy as a difference 
between the energy of the actual configuration and a configuration with very large and 
separated conductors. 



We find 




(10) 




(11) 
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3.1. The classical eigenvalues 

The solution of the Helmholtz equation in the annular region between eccentric cylinders 
has been considered in the framework of classical electrodynamics, fluid dynamics, and 
reactor physics, among others [201 EE] The eigenvalues have been computed using 
different approaches, as for instance conformal transformations that map the eccentric 
annulus onto a concentric one. As the two dimensional Helmholtz equation is not 
conformally invariant, the transformed equation has coordinate dependent coefficients 
and has to be solved numerically [22]- As is well known, it is very difficult to compute 
the Casimir energy from the numerical eigenfrequencies. It is more efficient to use 
the procedure outlined in Section 2, that only needs a function F with zeroes at the 
eigenvalues. Although for the eccentric annulus this function has been previously found 
in the literature [201 ETJ , for the benefit of the reader we include here a derivation of 
this result. 

The electromagnetic field inside an eccentric waveguide can be described in terms 
of TM and TE modes. The TM modes are characterized by a vanishing z component of 
the magnetic field, B z = 0. The other components of the electromagnetic field can be 
derived from the z component of the electric field, E z (r, 9, z, t) = E(r, Q^ e - lut + lzk ^ with 

E(r, 9) = £ [A m J m (Ar) + B m N m {\r)} e imd , (12) 

m 

where w 2 = k 2 z + A 2 , and (r, 9) are polar coordinates with origin at the center of the 
outer cylinder (see Fig. 1). One can also describe the z component of the electric field 
using polar coordinates (p, ip) with origin at the center of the inner cylinder (see Fig. 

1), 

E(p, p) = J2 [AA(Ap) + B„N n (\p)} e*"* . (13) 

n 

The perfect conductor boundary conditions imply that the z component of the electric 
field must vanish on the cylindrical shells 

A m J m (\b) + B m N m (\b)=0, 

A n J n {\a) + B n N n {\a) =0, (14) 

i.e., the functions E and E satisfy Dirichlet boundary conditions on the surfaces. The 
coefficients of the series in Eqs. (ll2jl and (jlHj) can be related to one another using the 
addition theorem for Bessel functions 

e^C m (Ap) = Y,e vpe C p (\r)J p ^ m (\e) , (15) 
v 

where C m denotes either J m or N m . Indeed, as at any point P in the annulus region one 
must have E(P) = E(P), it is possible to show that 

m 

B n = £5 m J n _ m (Ae). (16) 
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J n _ m (Ae) = 0. (17) 



Combining Eqs. (J14j) and (|16j) one obtains the linear, homogeneous system of equations 

Jn(Aa) _ J m (Xb) 
N n (Xa) N m (Xb)_ 

The solution of this linear system of equations is non trivial only if det[Q™] = 0, where 

Q™(a, b, e) = [J n (Xa)N m (Xb) - J m (Xb)N n (Xa)] J„_ w (Ae) . (18) 

This equation defines the allowed values for A, and therefore defines the eigenfrequencies 
of the TM modes. 

The TE modes can be treated in the same fashion. For these modes the z component 
of the electric field vanishes in the annulus region, E z = 0. The perfect conductor 
boundary conditions imply that the normal component of the magnetic field should 
vanish on the conducting shells, so now we must impose Neumann boundary conditions 
on the surfaces. The eigenvalues A for the TE modes are the solutions of det[Q™] = 0, 
where 

Ql E n (a, b, e) = [J' n (Xa)N' m (Xb) - J' m {Xb)N' n {Xa)} J n _ m (Ae). (19) 

In the concentric limit e = we have J n _ m (0) = 5 nm , the two matrices and 
Q™ become diagonal, and the equations for the eigenvalues are those of the concentric 
case jH] • In what follows we will use these matrices to define the function M that enters 
in Eq.(HH). 



3.2. The function M 

Roughly speaking, the function M that determines the Casimir energy through Eq. fTTf 
is the ratio of the function associated to the actual geometric configuration and the one 
associated to a configuration in which the conducting surfaces are very far away from 
each other. As the last configuration is not univocally defined, we will use this freedom 
to choose a particular one that simplifies the calculation. It turns to be convenient to 
subtract a configuration of two cylinders with very large (and very different) radii, but 
with the same eccentricity as that of the configuration of interest. 

We start considering the Dirichlet modes. To compute F lcy i(a) in Eq.® we note 
that the eigenfrequencies A for the geometry of a single cylinder of radius a surrounded 
by a larger one of radius R are defined by the equations 

J n (Xa) = , 

J n {Xa)N n {XR) - J n (XR)N n (Xa) = . (20) 

The first equation defines the eigenfrequencies in the region r < a and the second one 
gives the eigenfrequencies of the modes in the region a < r < R. Fi cy \(a) is the product 
of these two relations for all values of n, evaluated on the imaginary axis (A = iy = i(3/a). 
Namely, 

^icyi(a) = UUXa)[J n (Xa)N n (XR) - J n (XR)N n (Xa)} 

n 

= J(a)det[Q™(a,R,0)] , (21) 
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where we have introduced the notation J (a) = Y\ n J n (\a) to simplify the formulas below. 
The function F lcyl (oo) has the same expression, but replacing a by with R 1 very large 
but smaller than R. Using the asymptotic expansion of the modified Bessel functions it 
is easy to prove that F lcyl (a)/F lcy i(oo) ~ 2f3I n (f3)K n (f3)Ri/ a. The functions F and F^ 
in Eq.(jnj) are given by 

F = J{a)det[Q™{a, b, e)]det[Q™(6, R, 0)] 



J (a) 

m 



det[g iM (a,6,e)]F lcyl (6), 



Foo = ^|4det[Q™( J R 1 ,i? 2 ,6)]F lcyl (oo) 



(22) 
(23) 



where R\ < R 2 < R. As we already stressed, in order to define F^ we consider 
a configuration of two eccentric cylinders of large radii R± < R 2 and with the same 
eccentricity e of the original configuration. Evaluating the determinant in Eq. ()22|) on 
the imaginary axis one obtains 



det[Q iM (a,6,e)] = det 



-In 

7T 



P-) [K n {(3)I m {al3) 



" n I n {{3)K m {af3)} 



(24) 



Using again the asymptotic expansions of the Bessel functions one gets 
det[Q™(R 1 ,R 2 ,e)\ oc a/ n _ m (/3e/a)^ 

2y/R 1 R 2 p ' ^ e ec L ua ti° ns above can be combined 

to obtain 



M™(/3) = det 



x det/I 



t (/3e/a 



I n {a(3) 



-1) 



m+n 



I n {l3)K m {aft) 
K n (/3)I m (a(3) 



a 



(25) 



where I~^(fte/a) denotes the inverse matrix of I n - m {(3eJ a) and a = b/a. Computing 
explicitly the determinant one can show that the factor I m (aft)/I n (aft) cancels out. 
Moreover, writing M as a single determinant we get 

M™(/?)=det[5 np -A™], (26) 

with 



ATM _ ( i\n Iri(P) y^/ 1 n to K m (a(3) 

~ 1 > KJB) 4r l ' Ua0) 



13- 
a 



nip 



^n(ft) m Im(a(3) 

The addition theorem for the modified Bessel functions, C m (u ± v 
implies that I~p(/3e/a) 



a 



(27) 



A read 



.4 



TM 
np 



I n ((3) K m 
K n ((3) ^ I m ( 



l) m p I m _ p ([3e/a). Finally, the elements of the matrix 

aft) 



aft) 



ft- 



L p—m 



ft- 



(28) 



where we omitted a global factor 
determinant. 



■ l) n+p because it does not contribute to the 
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Figure 3. Exact Casimir interaction energy difference \AE\ between the eccentric and 
concentric configurations as a function of a = b/a for different values of 5 — e/a. Here 
AE = E\2 — E±2- Energies are measured in units of L/4na 2 . These results interpolate 
between the (a — 1)~ 5 behavior for a — ► 1, and the (a 4 log a) -1 behavior for a>l. 



The analysis for the TE modes is straightforward, the main difference being that 
Q TE (a, b, e) contains derivatives of those Bessel functions that do not depend on the 
eccentricity. Therefore, following similar steps it is possible to show that 



where 



M TE (/3)=det[<W-A™], ( 29 ) 

np ~ K'M \ I' m (aP) In ~ m \! j a) ^ \! j a) " (3 ° } 

The function M for the electromagnetic field is the product M = M TE M™, and 
therefore the interaction energy is the sum of the TE and TM contributions 

E 12 = — i \ d(3(3 InM(P) — - — ■= \ d(3 (3 lnM TE (/3) 
Ana 2 Jo Ana 2 Jo 

+ J— rd{3p\nM™([3) = E TE + E™. (31) 



Ana 2 

In order to calculate the exact Casimir interaction energy one needs to perform a 
numerical evaluation of the determinants. We find that as a approaches smaller values, 
larger matrices are needed for ensuring convergence. Moreover, for increasing values 
of the eccentricity e it is necessary to include more terms in the series defining the 
coefficients A™'™. In Fig. 3 we plot the interaction energy difference AE = E 12 — E^ 
between the eccentric (Ei 2 ) and the concentric (E^) configurations as a function of a for 
different values of 5. As we will show below, these numerical results interpolate between 
the PFA and the asymptotic behavior for large a. Fig. 4 shows the complementary 
information, with the Casimir energy as a function of S for various values of a, showing 
explicitly that the concentric equilibrium position is unstable. 
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Figure 4. Exact Casimir interaction energy difference AE between the eccentric 
and concentric configurations as a function of 8 = e/a for different values of a = b/a. 
Energies are measured in units of L/Aita 2 . The maximum at 6 = shows the instability 
of the concentric equilibrium position. 



3.3. Concentric Cylinders 



The exact Casimir interaction between concentric cylinders [6J can be obtained as a 
particular case of the exact formulas (J2%|) . (JHUjl . In the concentric limit e = 0, the 
matrices that appear in the definition of M TE and M™ become diagonal and the 
Casimir energy reads |Bj 

1 /3 ° d(3 p\nM cc (p), 



E 



12 



Ana 2 



where 



m cc (/5) = n 



1 - 



I n {(3)Kn{a(5) 



I' n {a(3)K>M 



(32) 



(33) 



(a(3)K n (/3)_ 

The first factor corresponds to Dirichlet (TM) modes and the second one to Neumann 
(TE) modes. 

The proximity limit a - 1 < 1 has already been analyzed for the concentric case 
jSj. In order to compute the concentric Casimir interaction energy in this limit it was 
necessary to perform the summation over all values of n. As expected, the resulting 
value is equal to the one obtained via the proximity approximation, namely 



E 



TE,cc 
12,PFA 



E 



TM.cc 
12.PFA 



-E 



EM,cc 
12,PFA 



7T 3 L 



1 



2""' rrJ1 720a 2 {a -If ^ 

Here i^ EM ' cc denotes the full electromagnetic Casimir interaction energy in the concentric 
configuration. 

In the large a limit one can show that only the n = term contributes to the 
interaction energy 



E cc as 



47T& 2 



dx x 



In 1 



jo(ggbfr) ' 
Jo(aO*o(i). 



+ ln 1 



.(35) 
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Figure 5. Modulus of the Casimir interaction energy in the concentric case as a 
function of a = b/a. Energies are measured in units of L/Aira 2 . These results 
interpolate between the (a — 1)~ 3 behavior for a — > 1, and the (a 2 log a) -1 behavior 
for a>l. 



Using the small argument behavior of the Bessel functions it is easy to prove that, in 
the limit a ^> 1, the TM mode contribution dominates, giving 

1.26L 



E 



12 



^_ rdxx^ 

4^6^ In a jo -^oW 



<n ■ (36) 
87ro z ma 

Note that the modulus of the energy decreases logarithmically with a. Fig. 5 depicts 
the exact Casimir interaction energy between concentric cylinders as a function of a, 
for values that interpolate between the above mentioned limiting cases. 

It is worth noticing that, while for small values of a both TM and TE modes 
contribute with the same weight to the interaction energy, the TM modes dominate in 
the large a limit. 



3.4- A cylinder in front of a plane 

It is interesting to see how the Casimir energy for the cylinder-plane configuration is 
contained as a particular case of the exact formula derived in Section 3.2. To do that 
let us consider a cylinder of radius a in front of an infinite plane, and let us denote by H 
the distance between the center of the cylinder and the plane. The eccentric cylinders 
formula should reproduce the cylinder-plane Casimir energy in the limit b, e — > oo 
keeping H = b — e fixed (see Fig. 1). 

We note that for x ^> h > 1 the ratio of Bessel functions appearing in Eq.()28|) can 
be approximated by 



lm - n \X) ^ 



m—p 



X) 



- p (x - h). 



(37) 



I m (x + h) 

This is trivially true for fixed m and large values of x, as can be seen from the large 
argument expansion of the Bessel functions. Moreover, using the uniform expansion of 
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the Bessel functions, it can be shown that Eq. (J37|) is also valid in the large m limit, for 
all values of x. Therefore we approximate 

^ ' — j 7 j rT^-^m— n(-^)-^m— — ^ ] K m (x ~\- /l)/ m _p_ n (s /i) 

= K n+p (2h), (38) 

where in the last equality we used the addition theorem of Bessel functions. Inserting 
this result (with x = (3e/a and h = (3H/a) in Eq.(j2HI) we get 

A™ - ^K n+p (2pH/a) = A™-"*, (39) 

which coincides with the known result for the Dirichlet matrix elements for the cylinder- 
plane geometry fT] ^] . The TE modes can be analyzed in the same fashion. Using 
that for large x 



K'Jx)I m (x) ^ 

K m {x)r m {x) 



(40) 



one can prove that 



X' {x + h) 

m 'mr ' ) m 



= -K n+p (2h). (41) 

Therefore, from Eq. ()30|) we obtain 

A™ - -^K n+p (2pH/a) = (42) 
which is the result for the TE modes in the cylinder-plane geometry [TT] I12j. 



4. Quasi-concentric cylinders 

We will now consider a situation in which the eccentricity of the configuration is much 
smaller than the radius of the inner cylinder, i.e., 8 — e/a <C 1. As discussed in jTl] 
this configuration may be relevant for performing null experiments to look for extra 
gravitational forces. Note that we do not assume that the radius of the inner and outer 
cylinder are similar, so the proximity approximation is in general not valid for this 
configuration. 

As described in the previous section, when e = the matrix defining the 
eigenfrequencies is diagonal. When considering a small non-vanishing eccentricity, the 
behavior of the Bessel functions for small arguments I m ^ n {(38) ~ {(35) n ~ m suggests 
that one only needs to use matrix elements near the diagonal. Using this idea, we 
will approximate the Casimir interaction energy by keeping only terms proportional to 
To, h and If. In this approximation, the matrices 5 np — A™ and 5 np — A^ p become 
tridiagonal matrices, and the e— dependent part of the Casimir energy will be quadratic 
in the eccentricity. 
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We will describe in detail the case of the Dirichlet modes; the treatment of Neumann 
modes is similar. To order 0(S 2 ) the only non- vanishing elements of the matrix A™ 
are 



A 



TM 



A 



TM 

n,n+l 



KM 



Kjap) r2/J . m , AT ra -i(q^) 2 | AT w+ i(a/3) 2 ' 



4TM _ 



'n+l 



l n {ap) l n -i 

K n (af3) | iT n+ i (<*/?) 

/„(«/?) J n+ i (a/5) 

(/?) |X(a/3) , K„ +1 (a/3) 



^n+l(/3) 



+ 



'n+l 



;«/?) 



'a/5) Jn+i 
I (6p)h(6P). 



a(3) 



(43) 



We split the matrix A™ into three terms, A™ = V™' cc + V™ +J\f™, where £>™> cc 
is the diagonal matrix corresponding to the concentric case, T>™ the diagonal part of 
the matrix that depends on 5, and A/"™ is the non-diagonal part of the matrix. The 
non- vanishing matrix elements are 

,TM,cc J n(P) K n (af3) 



•)TM 

n,n 

• /v n,n+l 



# n (/3) I n {a(3) 

TM -t~)TM,cc 



4 TM _ jy 

n,n n.n 



A™ 

■^n,n+l5 



TV™ 

- /v n+l,n 



^+V (44) 

Note that although T>™ = 0(5 2 ) and A/"™ = 0(5), both give quadratic contributions 
to the determinant. Up to this order we have 

-pTM 



In det [1 - A 



TMi 



In det [1 -D 1M ' CC ] + In det 
AT™ 



pTM,. 



+ In det 



1 - 



(45) 



I _ -pTM,cc 

The first term is associated to the interaction energy between concentric cylinders E 
studied in Section 3.3, and being ^-independent does not contribute to the force between 
eccentric cylinders. The second term can be easily evaluated 



•cc 
12- 



In det 



1 - 



V 



TM 



1 _ pTM,i 



In 1 



tr 



V 



TM 



1 _ pTM,i 



E 



v 



TM 



^-jTMjee 



.(46) 



To compute the last term in Eq. (j45|) we use that the determinant of an arbitrary 
tridiagonal matrix T of dimension p can be calculated using the recursive relation 
det[T{ p }] = T P)P det[T{p_i}] — T V)V -{T V ^ V det[T{ p _ 2 }], where T^} denotes the submatrix 
formed by the first k rows and columns of T. Up to quadratic order in S we obtain 
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(47) 
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Putting all together, the TM part of the Casimir interaction energy between quasi 
concentric cylinders can be written as 



77TM 
^12 



E 



TM,cc 
12 



Le 2 



4vra 4 



E/ 

„ Jo 



-, ^-.TM.cc 



^ n 



TM 



1-V 



TM,cc 
n+l,n+l 



.(4* 



Exact zero-point interaction energy between cylinders 



14 



Here 
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2 ' AK n {(3) 

I n ((3)I n+1 ((3) 
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(49) 



J„(a/3) J n+ i(a/3) 

The corresponding formulas for the TE modes can be obtained from these ones by 
replacing the Bessel functions by their derivatives with respect to the argument. 

The expression for the Casimir energy for quasi-concentric cylinders derived in 
this section is far simpler than the exact formulas Eqs. (j28|) . ([3U|h It is very useful for 
the analytical and numerical evaluation of the Casimir energy in the different limiting 
cases we will study below: the large distance limit (a <C b), for which one obtains 
a logarithmic decay of the energy, and the small distance limit (a ~ b), where the 
proximity approximation holds. The first case is very simple to handle because the 
energy is dominated by the lowest modes, while the second case is much more involved. 



4-1- Large distances: logarithmic decay 

When the ratio of the outer and the inner radii a = b/a is much larger than one, 
the exact Casimir energy is dominated by the lowest term n = in the summation. 
Moreover, it can be shown that the contribution of the Dirichlet modes is much larger 
than the contribution of the Neumann modes. Therefore, from Eq. ()31|) we get, in the 
limit a — > oo, 
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Using the small argument expansion of the Bessel functions it is easy to see that 
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(51) 



(52) 



I (x) 2a 2 \I (x) ■ h(x) 

In this expression, valid when a, e b, we kept the leading terms proportional to 
(a 2 In a) -1 and only the subleading terms that depend on the eccentricity. Inserting 
Eq.(|52|) into Eq.()50|) and computing numerically the integrals we find 



E 



12 



1.26 



SnbHnaV- 3 - 33 *)' (53) 

where the first term is the concentric contribution E^ ,cc derived before (see Eq.([36|)). 
It is worth to note that Eqs. and (p)Hj) have been derived under the assumption 

In a 3> 1, and therefore are valid for extremely large values of a. For intermediate values 
a ^> l,lna = 0(1), the interaction energy is also dominated by the Dirichlet n = 
term. The final result is still of the form given in Eq.([53|h with numerical coefficients 
that depend logarithmically on a. In Fig. 6 we plot the ratio between the exact 
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Figure 6. Ratio of the exact AE and asymptotic AEoo Casimir energy differences in 
the limit of small eccentricity e <C a. In the a — > oo limit, the Casimir energy difference 
between eccentric and concentric configurations decays logarithmically as (a 4 logo:) -1 . 



Casimir interaction energy difference AE = Ei 2 — E\\ and its asymptotic expression 
AEoo = E^ — E^ ,cc as a function of a. As mentioned, extremely large values of a are 
needed in order for the ratio of energies to asymptotically approach 1. From Eq. fjB^j) we 
see that the force between cylinders in the limit a, e <C b is proportional to Le/b A ln(b/a). 
The weak logarithmic dependence on the ratio b/a is characteristic of the cylindrical 
geometry (see also (H1H2]), an< ^ ^ * s a ^ so f° un d in the electrostatic counterpart of the 
Casimir energy, that we briefly analyze next. 

The electrostatic capacity for the system of two eccentric cylinders is given by 

C = 2 ™f , (54) 

where Y = (a 2 + b 2 — e 2 )/2ab, and eo is the permittivity of vacuum. Therefore, the 
electrostatic force between cylinders kept at a fixed potential difference V is 

F _ ± 7re V 2 L 

^ elec " ab VY^T \n 2 [Y + VY^l] ' 1 j 

In the quasi-concentric case we can set e = in the definition of Y. In the large a limit 
we get 

Le 



elec 



(56) 



Just as in the Casimir case, in the limit a<5 the Coulomb force decays logarithmically 
with the ratio a/b. 

4-2. Small distances: the proximity approximation 

The proximity limit for concentric cylinders has been reviewed in Sec. 3.3; the case of 
a cylinder in front of a plane has been considered in detail in In this section we 
extend these results to the case of quasi-concentric cylinders. 




Figure 7. Ratio of the exact and PFA Casimir interaction energy differences 
AE = E12 — between eccentric (-E12) and concentric cylinders in the limit 

of small eccentricity e <C a. The curve EM denotes the full electromagnetic Casimir 
energy. 



We will concentrate on calculating the Casimir interaction energy difference AE = 
E\2 — E^2 between the eccentric (-E12) and the concentric (-E12) configurations. As the 
small distance limit is dominated by the large-n modes, the key point in the derivation 
of the PFA from the exact expression of the Casimir energy is the use of the uniform 
approximation for the Bessel functions. In the large n limit, and to leading order in 
a — 1 one has 

I n (f3) K n {a(3) 



K n {(5) I n (a(3) 
IM K n±l {a(3) 
K n ((3) I n±1 (a(3) 



-2n(a-l)h(x) 



-2n(a-l)h(x) 



1 + h(x) 



X 



±2 



(57) 



where (3 = nx and h(x) = yl + x 2 . Inserting these approximations in Eq. (j49jl we get 
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The contribution to the interaction energy coming from the diagonal part of the 
matrix can be written as (see Eq. 

LS 2 
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(59) 



n,k>l 



where we replaced the sum over all integers n by twice the sum over the positive integers 
(the term n = gives a sub leading contribution for small a — 1). Inserting the uniform 
expansions into Eq. (j59j) and changing variables in the integral we obtain 

L5 2 



AE™ 



4na 2 
1 + 



E 



n 



dx x 3 e~ 2n ( a ~^ h ^ k 



n,k>l 

1 fl + h(x) 

2 1 x 



+ 



1/1 + h(x) 



-2 



x 



(60) 
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To leading order in a — 1 the sum over n gives J2n n ^ e 2n<yOL 1 ' h ^ x ' — 24/[2/t(x)(a — 1)£;] 5 . 
Next we perform first the sum over k and then the integral over x. We get 

3 L5 2 ((5) 



AE™ 



i\a 2 {ct — l) 5 ' 



(61) 



where ((x) is the Riemann zeta function. 

The evaluation of the non diagonal contribution to the Casimir energy can be done 
using similar steps, starting from Eq. (}4Tj) . In the proximity limit, we can approximate 

lfn+i by ^™' €C m the denominator of that equation. Therefore, using Eq.(0HJ)we 
write the non diagonal contribution to the energy as 

- ~ dp P 3 E AC™ k {Vlf-r 1 - (62) 
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Now we use the uniform expansion for the Bessel functions in Eq. (J49|) to obtain 
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(63) 



Replacing Eq. (j63|) into Eq. (j62|) we get 
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(64) 



As before, we first compute the sum over n, and expand the result to leading order in 
a — 1. The sum over k can be calculated using that J2k>i (k+i) 5 = C(4) — C(5). Finally, 
we compute analytically the remaining integrals to get 

3 L5 2 
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The contribution of the Dirichlet modes to the Casimir interaction energy in the limit 
a — > 1 is therefore 



AE 



TM 



AE™ + AE 



'TM 
'MV 



L6 2 



7T" 



(66) 



Jv ' ~~ N ' U a 2 (a- 1) 5 240' 
It can be shown that the contribution of the TE modes to the interaction energy in the 
short distance limit is equal to that of the TM modes, as expected from the parallel 
plate configuration. Indeed, the uniform expansion for the ratio of Bessel functions is 
equal to the expansion for the derivatives, i.e., 
I'M K' n (aP) 
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and therefore all calculations can be repeated without changes. 

The final result for the Casimir interaction energy difference in the small distance 
approximation is 

1 7T 3 Le 2 
AE™ A = AE™ A = -A£™ = - 240fl4(a _ 1)5 , (68) 

where AE EM denotes the full electromagnetic Casimir energy difference between 
eccentric and concentric configurations. Fig. 7 depicts the ratio of the exact 
Casimir energy difference AE and the PFA limit for the almost concentric cylinders 
configuration. As evident from the figure, PFA agrees with the exact result at a few 
percent level only for a close to unity, and then it noticeably departs from the PFA 
prediction. The resulting PFA expression for the attractive Casimir force between quasi- 
concentric cylinders reads 



1>KY_« 3 - (6g) 



n" ( 1. 



60a 4 (a-l) 5 ' 
that reproduces the result previously obtained in [TJ] 



5. Conclusions 



We have derived an exact formula for the Casimir interaction energy between eccentric 
cylinders using a mode summation technique. This formula is written as an integral 
of the determinant of an infinite dimensional matrix, and it reproduces as a particular 
case the interaction energy between concentric cylinders, and as a limiting case the 
energy in the cylinder-plane geometry. In the quasi-concentric case, the infinite 
dimensional matrix becomes tridiagonal, and hence much easier to deal with than the 
exact formula when performing analytic and numerical calculations. We have carried 
out the numerical evaluation of the Casimir interaction energy using both the exact and 
tridiagonal formulas, and studied different limiting cases of relevance for Casimir force 
measurements. 

The large and small distance limits were analyzed. In the former case, the Casimir 
energy is dominated by the lowest modes, and shows a weak logarithmic decay, typical of 
cylindrical geometries. In the latter case, the Casimir energy is dominated by the highest 
modes, and the exact formula reproduces the proximity approximation. We found that 
the first order correction (a — 1 <C 1) to PFA for the quasi-concentric cylinders has 
the form AE / AE^-pA = 1 + s(ct — 1) + 0((a — l) 2 ), where the coefficient of the linear 
curvature correction is positive, s > 0, both for TE and TM modes. This contrasts with 
the first order corrections to PFA in the cylinder-plane configuration, where the linear 
curvature correction to TM modes is positive, while the one for TE modes is negative 
P. 

The exact Casimir force computed in this paper, in particular for the quasi- 
concentric configuration, offers a qualitatively different approach for implementing new 
experiments to measure the Casimir force and to search for extra-gravitational forces 



Exact zero-point interaction energy between cylinders 



19 



in the micrometer and nanometer scales, since it opens the possibility of measuring the 
derivative of the force using (Cavendish-like) null experiments. 
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